On the 1/ / Frequency Noise in Ultra-Stable 
Quartz Oscillators 

Enrico Rubiola 3 , and Vincent Giordano 

3 web page http://rubiola.org 



sciences & tecf^logies 



FEMTO-ST Institute 
CNRS and Universite de Franche Comte, Besangon, France 



2nd February 2008 



Abstract 

The frequency flicker of an oscillator, which appears as a l// 3 line in 
the phase noise spectral density, and as a floor on the Allan variance plot, 
originates from two basic phenomena, namely: (1) the 1/f phase noise 
turned into 1/f frequency noise via the Leeson effect, and (2) the 1/f 
fluctuation of the resonator natural frequency. The discussion on which is 
the dominant effect, thus on how to improve the stability of the oscillator, 
has been going on for years without giving a clear answer. This article 
tackles the question by analyzing the phase noise spectrum of several 
commercial oscillators and laboratory prototypes, and demonstrates that 
the fluctuation of the resonator natural frequency is the dominant effect. 
The investigation method starts from reverse engineering the oscillator 
phase noise in order to show that if the Leeson effect was dominant, the 
resonator merit factor Q would be too low as compared to the available 
technology. 
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Symbol list 



A amplifier (voltage) gain 

bi coefficient of the power-law representation of S v (/) 

b(t) resonator impulse response in the phase space 

B(s) Laplace transform of b(t) 

f Fourier frequency (near dc) 

f c amplifier corner frequency, divides white and 1// regions 

Jl Leeson frequency 

F amplifier noise figure 

hi coefficient of the power- law representation of S y (f) 

H(s) loop transfer function (Laplace transform) 

H(s) loop transfer function (Laplace transf.) in the phase space 

j imaginary unit, j = —1 

k 1.38x10—23 J/K, Boltzmann constant 

P, Po carrier power. Also P , Pb, P m , etc. 

Q resonator merit factor 

R frequency-srability ratio 

s complex frequency s = a + ju, in Laplace transforms 
S(f), S x (f) single-sided power spectrum density (of the quantity x) 

t time 

T, T absolute temperature, reference temperature (T = 290 K) 

v(t) (voltage) signal, as a funtion of time 

Vo, Vi peak carrier voltage 

y fractional frequency fluctuation, y = 

a(t) fractional amplitude fluctuation 

0(s) resonator impulse response (Laplace transform) 

v, v a frequency, carrier frequency 

Cj/(r m ) Allan deviation of the quantity y 

r m as in <J y {T m ), measurement time 

t resonator relaxation time 

<p(i) oscillator phase fluctuation 

$(s) Laplace transform of ip(t) 

tp(t) amplifier phase fluctuation 

'I'(s) Laplace transform of ip(t) 

u) angular frequency 



Note: u> is used as a shorthand for 2ttu or 2nf, and viceversa 
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1 Introduction and summary 

In the domain of ultra-stable quartz oscillators used in the most demanding ap- 
plications, like space and atomic fountain clocks, we notice that the frequency 
flicker is often the most critical parameter. The required stability is some- 
times in the upper 10 -14 (Allan deviation) at 1-30 s or so, which can only be 
achieved in the lower HF band (5-10 MHz), and after selection. In such cases, 
identifying the dominant flicker mechanism is far from being trivial. Whereas 
some authors strongly suggest that the amplifier noise can be the parameter 
that limit the frequency stability, rather than the flickering of the resonator 
natural frequency |Wal95l lBMG + 99~] , the general literature seems not to give 
a clear answer. This conclusion results from a set of selected articles, which 
includes the measurement of the frequency stability WW75 RGBG00 and the 
interpretation of the flicker noise of crystal resonators |Kro 98 Kro05j; the de- 
sign fundamentals of the nowadays BVA resonators |Bes77) : some pioneering 
works on the low-frequency noise in quartz oscillators |BOM75| IDri75j ; more 
recent articles focusing on specific design solutions for ultra-stable oscillators 
|Nor91l INor96l ICCG+981 ICCL+031 IT JA97| : and, as a complement, a thorough 
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review of the SiC-2 crystal for the resonator fabrication is found in [Bri85 . Con- 
versely, in everyday-life oscillators, which span from the low-cost XOs to the 
OCXOs used in telecommunications and instrumentation, the relative simplic- 
ity of the low-noise electronics required indicates that the frequency flicker is 
chiefly the 1/ f fluctuation of the resonator. 

In a previous work |Rub05j . now extended to more commercial products 
and laboratory prototypes, we have analyzed the phase noise spectrum of some 
oscillators, aiming at understanding the internal mechanisms and parameters. 
We look at the phase-noise spectrum from the right hand to the left, hence 
from the higher Fourier frequencies to the lower, matching theory, technology 
and physical insight. In this way we get information on the sustaining amplifier 
on the output buffer, on the Leeson effect and on the resonator. 

In this article we first explain the phase noise mechanisms in amplifiers. 
Then we introduce the Leeson effect, which consists of the phase-to-frequency 
conversion of noise below the resonator cutoff (Leeson) frequency Jl — ^q- Fi- 
nally, we analyze the phase noise spectral density S ip (f) of a few oscillators. 
The conclusion that the resonator natural frequency is the main cause of fre- 
quency flickering is based on experimental facts. After taking away the effect 
of the output buffer, we calculate the frequency f'£ at which the oscillator f~ 3 
line crosses the / _1 line of the sustaining amplifier. Provisionally assuming 
that f'l is the the Leeson frequency, we observe that the resonator merit factor 
Q s = ^jtt thereby calculated is far too low for a high-tech resonator. Conversely, 
under any reasonable assumption about the true merit factor, the Leeson effect 
is found at a frequency /l <§C f£. Therefore the Leeson /~ 3 line on the S v (f) 
plot is well hidden below the resonator fluctuation. 



2 Phase noise fundamentals 

Let the quasi-perfect oscillator sinusoidal signal of frequency i/q 

v(t) = V [l + a(i)] cos[27w i + <p(t)] . (1) 

where a(t) is the fractional amplitude noise, and Lp(i) is the phase noise. The 
AM noise is not essential to this work. The phase noise is best described in 
terms of S v (f), i.e., the one-sided power spectral density of <p(t), as a function 
of the Fourier frequency /. In addition to /, we use the angular frequency lo for 
both carrier-related frequencies (us = 2ttv), and Fourier frequencies {lo — 2tt/) 
without need of introducing it, and the normalized frequency fluctuation y = 
v ~" a . The quantities v, f and y refer to one-sided transforms, to to two-sided 
transforms. Frequency fluctuations are described in terms of S y (f), related to 



s y (f) = ^s v (f) 



(2) 
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C4 




i f depends on P 

Figure 1: Typical phase noise spectral density of an amplifier. 



A model that has been found useful in describing the oscillator noise spectra 
is the power-law 

2 

Sy(f) = E W S M) = E b ^ ■ ( 3 ) 

i=-2 i=— 4 

Our main concern is the frequency flickering term 6-3/~ 3 , which is related to 
the Allan variance by 

&y = 21n(2)ft_! = 21n(2) ^ , (4) 

constant, i.e., independent of the measurement time. 

Finally, the general background on phase noise and frequency stability is 
available from numerous references, among which we prefer Rut 78 , CCI90 , 
|Kro83| , and jVA89l Vol.1, Chapter 2]. A IEEE standard is also available 
|Vi99| . 

3 Phase noise in rf (and microwave) amplifiers 
3.1 White noise 

The equivalent noise spectrum density at the amplifier input is N = FkTo, 
where F is the noise figure and kTo is the thermal energy. This type of noise 
is additive. In the presence of a carrier of power Pq, the phase noise spectral 
density is 



S v (f) = b o.f° (constant) 



(5) 
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with 



When amplifiers are cascaded, the noise contribution of each stage is divided by 
the gain of all the preceding stages (Friis formula |Fri44| ). Accordingly, in most 
practical cases the total noise is chiefly the noise of the first stage. Of course, 
this also holds for phase noise. 



3.2 Flicker noise 

Understanding the close-in noise starts from the bare observation that the out- 
put spectrum is of the white type — flat in a wide frequency range — when the 
carrier power is zero, and that noise shows up close to the carrier only when a 
sufficiently large carrier power is present at the amplifier output. The obvious 
consequence is that the close-in flickering results from a parametric effect by 
which the near-dc flicker noise modulates the carrier in amplitude and phase. 

The simplest way to understand the noise up-conversion is to model the 
amplifier signal as a nonlinear function truncated to the 2nd order 

v (t) = aiv t (t) + a 2 v 2 t {t) + . . . , (7) 

in which the complex input signal 

v i (t) = V i e^° t +n'(t)+jn"(t) (8) 

contains the carrier and the internally generated near-dc noise n(t) = n'it) + 
jn"(t). Rather than being an easy-to-identify voltage or current, n(t) is an 
abstract random signal that also accounts for the efficiency of the modulation 
process. Combining Q and JHJ and selecting the terms close to the carrier 
frequency uiq, we get 

v (t) = V l {a 1 + 2a 2 [n'(t)+jn"(t)]}e^ at . (9) 

Hence, the random fluctuations are 

a{t) = 2 — n'{t) and ip(t) = 2 — n"(t) . (10) 

Deriving Eq. {XSjl , the statistical properties oin'(t) and n"(t) are not affected 
by the carrier power. This accounts for the experimental observation that the 
amplifier phase noise given in rad 2 /Hz is about independent of power in a wide 
range | HWBri8l IWF.lHfl IHHWW 03 . Thus 

S v (f) = b-J- 1 6_i re constant . (11) 

Of course, some dependence on Pq remains. We ascribe it to terms of order 
higher than 2 in Q , and to the effect of the large signal regime on the dc bias. 
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In the case of bipolar amplifiers used in HF/VHF amplifiers, 6_j is in the range 
of I(T 12 to fCT 14 rad 2 /Hz (-120 to -140 dBrad 2 /Hz). 

When m amplifiers are cascaded, the The Friis formula does not apply. 
Instead, the phase noise barely adds 

m 

(&-l)cascadc = ^(6-l)i . (12) 

i=l 

This occurs because the 1// phase noise is about independent of power. Of 
course, the amplifiers are supposed independent. 

3.3 Phase noise spectrum 

Combining white noise [Eq. (JjJJ] and flicker noise [Eq. (HJ], there results the 
spectral density S v (f) shown in Fig.^ It is important to understand that the 
white noise term 60 / depends on the carrier power Pq, while the flicker term 
does not. Accordingly, the corner frequency f c at which = b is 

a function of P , thus f c should not be used to describe noise. The parameters 
F, and P should be used instead. 

4 Phase noise in feedback oscillators 

4.1 The Leeson effect 

Figure|2]shows a model for the feedback oscillator, and its equivalent in the phase 
space. All signals are the Laplace transform of the time-domain quantities, as a 
function of the complex frequency s — a + juj. The oscillator transfer function 
is derived from Fig. 0A according to the basic rules of linear systems 

"w^zh™ (13) 

A0(') 

Stationary oscillation occurs at the angular frequency loq at which Af3(Ju>) = 1, 
thus \A/3(juj)\ = 1 and &rg[Af3(jui)) = 0. This is known as the Barkhausen 
condition for oscillation. At s — jujQ the denominator of H(s) is zero, hence 
oscillation is sustained with zero input signal. Oscillation starts from noise or 
from the switch-on transient if ^{A/3(s)\ s —j Uo } > 1 (yet only slightly greater 
than 1 for practical reasons). When the oscillation reaches a threshold ampli- 
tude, the loop gain is reduced to 1 by saturation. The excess power is pushed 
into harmonics multiple of u>q, and blocked by the resonator. For this reason, 
at luq the oscillator operates in quasi-linear regime. 

In most quartz oscillators, the sustaining amplifier takes the form of a nega- 
tive resistance that compensates for the resonator loss. Such negative resistance 
is interpreted (and implemented) as a transconductance amplifier that senses 
the voltage across the input and feeds a current back to it. Therefore, the 
negative-resistance oscillator loop is fully equivalent to that shown in Fig. [21 
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P(s) 

resonator 



oscillator buffer 



B - Phase-space equivalent 

>P(s) <J>(s) 



oscill 
out 



B(s) 



*0« 



mam 
out 



oscillator buffer 



Figure 2: Oscillator model and its phase-space equivalent. For the sake of 
simplicity all the dependence on s is moved to /3(s) , hence the gain A is assumed 
constant. The scheme emphasizes the amplifier phase noise. Amplitude noise is 
not considered. 



In 1966, D. B. Leeson |Lee66| suggested that the oscillator phase noise is 
described by 



,2 1 



S^,(f) (Leeson) 



(14) 



This formula calls for the phase-space representation of Fig.[2]B, which deserves 
the following comments. 

The Laplace transform of the phase of a sinusoid is probably the most com- 
mon mathematical tool in the domain of PLLs |KF72I IHar79l IBes99l |Ega98| . 
Yet it is unusual in the analysis of oscillators. The phase-space representation 
is interesting in that the phase noise turns into additive noise, and the sys- 
tem becomes linear. The noise-free amplifier barely repeats the input phase, 
for it shows a gain exactly equal to one, with no error. The resonator transfer 
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function, i.e., the Laplace transform of the impulse response, is 

B(s) = - J— t = ^ . (15) 

The inverse time constant is the low-pass cutoff angular frequency of the res- 
onator 

^ = Z = 2Q ■ 



r 



The corresponding frequency 

f = ^ _ — = 

71 2tt 27tt 2Q V ' 

is known as the Leeson frequency. Equation (|15[) is proved in two steps: 

1. Feed a Heaviside step function nU(t) in the argument of the resonator 
input sinusoid. The latter becomes cos [u>ot + KU(t)]. 

2. Linearize the system for k — ► 0. This is correct in low phase noise condi- 
tions, which is certainly our case. Accordingly, the input signal becomes 
cos(wo£) — K sm(Lu t)U (t) . 

3. Calculate the Laplace transform of the step response, and use the prop- 
erty that the Laplace transform maps the time-domain derivative into a 
multiplication by the complex frequency s. The Dirac function S(t) is the 
derivative of U(t). 

The full mathematical details of the proof are available in |Rub05l Chapter 3] . 

Applying the basic rules of linear systems to Fig. [2JB, we find the transfer 
function 

n{s) = W) = t^W) - ' (18) 



thus 



(19) 



The Leeson formula (|14J) derives from Eq. (|19|l by replacing 

Lu = 2nf and r = — . (20) 

The transfer function Ti{s) has a pole in the origin (pure integrator), which 
explains the Leeson effect, i.e., the phase-to-frequency noise conversion at low 
Fourier frequencies. At high Fourier frequencies it holds that Ti(juS) = 1. In 
this region, the oscillator noise is barely the noise of the sustaining amplifier. 
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A - High Q, low v (xtal) B - Low Q, high v Q (microw.) 




Figure 3: Oscillator phase noise spectra, not accounting for the output buffer. 

The amplifier phase noise spectrum contains flicker and white noise, i.e., 
s i>(f) = ( & -i)am P ii f^ 1 + (M amp ii- Feeding such S^(f) into the Leeson formula 
H14|l. the oscillator S v> (f) can only be one of those shown in Fig. [31 Denoting with 
f c the corner frequency at which flicker noise equals white noise, we often find 
/l < fc in HF/VHF high-Q oscillators, and fa > f c in microwave oscillators. 
In ultra-stable HF quartz oscillators (5-10 MHz), the spectrum is always of the 
type A (f L < f c ). 

4.2 Output buffer 

The phase noise S-^b(f) of the output buffer barely adds to the oscillator phase 
noise 

S(po\f) = 

This a consequence of the flicker noise mechanism explained Section 13.21 and 
inherent in the model of Fig.[3B. 

4.3 Resonator stability 

The oscillator frequency follows the random fluctuation of the resonator natural 
frequency. However complex or tedious the formal proof for this statement can 
be, the experimentalist is familiar with the fact that the quartz oscillator can 
be frequency-modulated by a signal of frequency far higher than the Leeson 
frequency. For example, a 5 MHz oscillator based on a Q = 2xl0 6 resonator 
shows a Leeson frequency of 1.25 Hz (see Table EJ, while it can be modulated 
by a signal in the kHz region. Additionally, as a matter of fact, the modulation 
index does not change law from below to beyond the Leeson frequency. This 
occurs because the modulation input acts on a varactor in series to the quartz, 
whose capacitance is a part of the motional parameters. 



1 + 

P 4Q 2 



S^(f) + S^ b (f) 



(21) 
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4.4 Other effects 

The sustaining amplifier of a quartz oscillator always includes some kind of feed- 
back; often the feedback is used to implement a negative resistance that makes 
the resonator oscillate by nulling its internal resistance. The input admittance 
Yi seen at the amplifier input can be represented as 

Y l = + y/ r) , (22) 

that is, the sum of a virtual term (v) plus a real term (r). The difference 
between 'virtual' and 'real' is that in the case of the virtual admittance the input 
current flows into the feedback path, while in the case of the real admittance 
the input current flows through a grounded dipole. This is exactly the same 
concept of virtual impedance routinely used in the domain of analog circuits 
|Fra97l Chapter 1]. The admittance Y^ also includes the the effect of the 
pulling capacitance in series to the resonator, and the stray capacitances of the 
electrical layout. As a consequence, the fluctuation 5Y^ is already accounted 
for in the amplifier noise, hence in the model of Fig. |2 while the fluctuation 
SY^ is not. On the other hand, K interacts with the resonator parameters, 
thus 5Y^ yields frequency fluctuations not included in the Leeson effect. The 
hard assumption is made in our analysis, that |J1^| 2 <C |£f/^| 2 . In words, 
we assume that the fluctuation of the electronics are chiefly due to the gain 
mechanism of the amplifier. Whereas the variety of circuits is such that we can 
not provide a proof for this hypothesis, common sense suggests that electronics 
works in this way. 



5 Analysis of the oscillator phase noise 

This section address the core question, wether the l// 3 noise observed on the 
oscillator S v (f) plot is due to the Leeson effect, or it originates in the resonator. 
The interpretation method is shown in Fig. 0] and discussed below. 

1. We start from the spectrum, measured or taken from the oscillator spec- 
ifications. The first step is to remove the residual of the mains (50 or 60 
Hz and multiples) and other stray signals, and to fit the spectrum with 
the power-law [Eq. @]. This process is called parametric estimation of 
the spectrum. With a pinch of experience, sliding old-fashion squares on 
a A4-size plot gives unexpectedly good results. Otherwise, the mathe- 
matical methods explained in |.TW68I IPW98) are useful. After this, the 
spectrum looks like that of Fig. 0]A, 

2. The term &o/° is chiefly due to the sustaining amplifier, hence the amplifier 
input power can be calculated using Eq. © 

P„ - . (23) 

In the absence of information, it is wise to take F = 1.26 (1 dB). To the 
extent of our analysis, estimating Pq is only a check of plausibility. 



Table 2: Estimated Parameters of some Ultra-Stable Oscillators. 



Oscillator 



^0 



(i>-3)tot (&-l)tot (&-l)amp /I /£' 



Qs 



Qt 



h (6-a)t 



Ref.& 
Note 



Oscilloquartz r 
8600 '* 



-124.0 



-131.0 



-137.0 2.24 4.5 5.6xl0 5 1.8xl0 6 1.4 -134.1 10.1 



TJA97ll5sT| 

(1) 



Oscilloquartz r 
8607 ' 



-128.5 



-132. 



-138.5 1.6 3.2 7.9xl0 5 2xl0 6 1.25 



-136.5 



|TJA971l5sTj 

(1) 



CCG+98 CC L+03llcnT5j 

(1) 



CMAC 
Pharao 



-132.0 



-135.5 



-141.1 1.5 3 8.4xl0 5 2xl0 6 1.25 -139.6 7.6 



FEMTO-ST 
LD protot. 



10 



-116.6 



-130.0 



-136.0 4.7 9.3 5.4xl0 5 1.15xl0 6 4.3 



-123.2 6.6 



GSJBM04 

(3) 



Agilent 
10811 



10 



-103.0 



-131.0 



-137.0 25 50 1x10 s 



7xl0 5 7.1 -119.9 16.9 l^ W81 1 



Agilent 
prototype 



10 



-102.0 



-126.0 



-132.0 16 32 1.6xl0 5 7xl0 5 7.1 



-114.9 12.9 

(5) 



TtTnfn^n 100 -67.0 -132? -138? 1800 3500 1.4xl0 4 8xl0 4 625 -79.1 15.1 
501-04623 (3) 

Z ^ dB dB dB ~ ~ ' " ' " ~ dB ~ 
unit MHz rad2/Hz rad2/Hz rad2/Hz Hz Hz (none) (none) Hz dB 

Notes 

(1) Data are from specifications, full options about low noise and high stability. 

(2) Measured by CMAC on a sample. CMAC confirmed that 2xl0 6 < Q < 2.2xl0 6 in actual conditions. See Fig. El 

(3) LD cut, built and measured in our laboratory, yet by a different team. All design parameters are known, hence Qt- 

(4) Measured by Hewlett Packard (now Agilent) on a sample. 

(5) Implements a bridge scheme for the degeneration of the amplifier noise. Same resonator of the Agilent 10811. 

(6) Data are from specifications. See Fig. HJ1 
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3. Feeding the oscillator 6_ 3 term into Eq. (@J, we calculate the floor of the 
Allan deviation <j y . We check on the consistency between calculated value 
and specifications or measures, if available. 

4. At first sight, the cutoff frequency f' L (Fig. 0JA) can be taken for the 
Leeson frequency because there the slope changes from /~ 3 to Yet 
the term 6_i/~ 1 contains the flicker of the sustaining amplifier and of the 
output buffer, which add [Equations (|12|l and 121(1 ], For this reason, f' L 
can not be the Leeson frequency. 

5. Actual oscillators have 2-4 buffer stages, the main purpose of which is 
to isolate the feedback loop from the environment in order to ensure fre- 
quency stability and to prevent injection locking. Owing to the Leeson 
effect, a wise designer will spend the lowest- noise technology in the sus- 
taining amplifier, rather than in the buffer. Thus, we assume that the 
buffer contributes 3/4 of the total noise, and that sustaining amplifier con- 
tributes 1/4 (—6 dB). Accordingly, we plot the line fe-ia/ 1 in Fig.Q]B, 
6 dB below the total flicker. 

6. After taking away the buffer noise, the continuation of the bsf~ 3 line 
meets the 6-i a / _1 line at f — f'£. The latter is a new candidate for the 
Leeson frequency. Feeding f'l into Eq. I(17|l . we calculate the resonator 
merit factor Q s (the subscript s stands for 'spectrum') 

Qs = Wl ■ (24) 

7. Technology suggests a merit factor Q t (the subscript t stands for 'technol- 
ogy') significantly larger than Q s , even in actual load conditions. Feeding 
Qt into Eq. (|17|) . we calculate /l based on the actual merit factor 

h 'W.' (25) 

as shown in Fig. Q]C. There follows a phase noise term (6-3)^, which 
account for the Leeson effect alone. 

8. Given Q t 3> Q s , thus Jl <C f'l, the Leeson effect is hidden. Consequently, 
the oscillator / _3 phase noise is chiefly due to the fluctuation of the res- 
onator natural frequency. 

We introduce the stability ratio R, defined as 

R = ( ^ )osci " (floor), (26) 

[P^y J Leeson 

and related to the other oscillator parameters by 



„ / (&-3)tot Qt f L 

R - vrar = oi = h' (7) 
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Figure 5: Phase noise of the CMAC Pharao 5 MHz quartz OCXO. Courtesy of 
CMAC. Interpretation and mistakes are of the authors. 



This can be demonstrated from the 6_3 term of the Leeson formula 1)14(1. using 
Equations J3J and JTJJ. The parameter R states how bad is the actual oscillator, 
as compared to the same oscillator governed only by the Leeson effect, with the 
resonator fluctuations removed. Thus, R = 1 (0 dB) indicates that the oscillator 
f~ 3 phase noise comes from the Leeson effect. Equal contribution of resonator 
and Leeson effect yield R = \J~2 (3 dB), while R 3> y/2 is found when resonator 
instability is the main cause of /~ 3 phase noise. In all cases we have analyzed, 
discussed in the next Section, we find R of the order of 10 dB, with a minimum 
of 6.6 dB. This means that the Leeson effect is hidden below the frequency 
fluctuation of the resonator. 

Coming back to the estimation of the 1// noise of the sustaining amplifier 
it is to be remarked that if the 1/ / noise of this is lower than 1/4 of the total 
flicker, f'£ is further pushed on the right hand on Fig. 0]B-C, which reinforces 
the conclusion that the resonator is the main cause of frequency fluctuation. 



6 Experimental data and discussion 

Figure shows the phase noise spectrum of a 5 MHz oscillator, out of a small 
series intended as the flywheel for the space Cesium fountain clock Pharao 
| LLS+01llBLA+04| . On this plot, the reader can follow the interpretation pro- 
cess explained in Section [SJ and illustrated in Fig.0] Guessing on technology, 
the merit factor was estimated to be 2xl0 6 . Afterwards, the manufacturer con- 
firmed |Canj that Qt is between 2xl0 6 and 2.2xl0 6 in actual load conditions 
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for that series of oscillators, and that the flicker noise of the sustaining amplifier 
is less than 1/4 (—6 dB) of the total flickering. This validates our conclusions. 

Table |21 shows the results of our analysis on some oscillators. The ability 
to estimate the resonator merit factor is necessary to understand the oscillator 
inside. Experience indicates that the product uqQ is a technical constant of the 
piezoelectric quartz resonator, in the range from lx 10 13 to 2x 10 13 . As a matter 
of fact, the highest values are found in the 5 MHz resonators. In load conditions, 
the resonator merit factor is somewhat lower. The actual value depends on 
frequency, on the designer skill, and on the budget for implementation. A bunch 
of data are available from |OB85l IKro05l Wal95 , and from our early attempts 
to measure the resonator frequency stability |RGBG00| . The oscillate rs we 
have analyzed exhibit the highest available stability, for we are confident about 
published data. The Agilent 10811 (hence the Agilent prototype) is closer to 
the routine production, and probably closer to the cost-performance tradeoff, 
as compared to the other ones, thus understanding oscillator the inside is more 
difficult. Nonetheless, in this case the value of Q s is so low that there is no 
doubt that it can not be the resonator merit factor. 

In the case of the Oscilloquartz 8607, the /~ 3 noise is too low for it to be 
extracted from the S v (f) spectrum available on data sheet, which starts from 
1 Hz. Yet, we can use the device specifications S v (f)\xn z = —127 dBrad 2 /Hz, 
SVC/)|iohz = -142 dBrad 2 /Hz, and S v (f)\ lkHx = -153 dBrad 2 /Hz. In fact, 
looking at the spectrum and at the Allan variance it is clear that at / = 1 Hz and 
/ = 10 Hz the terms &_3/" 3 and determine S^lf), with at most a minor 

contribution of bo. It is also clear that 5 ¥> (/)|ikHz — &o- Thus 6_3 and 6_i are 
obtained by solving a system of two equations like S v (f) = b-3f~ 3 +b-if~ 1 +bo, 
at 1 Hz and 10 Hz. 

In the case of the Wenzel 501-04623 oscillator (Fig. |5J), the specifications 
available on the manufacturer web site consist of a few points, while the whole 
spectrum is not published. Experience indicates that in the case of 100 MHz 
oscillators the / _1 line tends to be hidden by the frequency flickering. That 
said, we can only guess that the / _1 noise of the sustaining amplifier is similar 
to that of other oscillators. This is sufficient to estimate /£', and to notice that 
the merit factor Q s is far too low as compared to the state of the technology, and 
to conclude that the f~ 3 phase noise is due to the fluctuation of the resonator 
natural frequency. It is to be remarked that the power at the amplifier input 
is of the order of 10-20 /xW in all other cases, and of 1 mW here. In addition, 
the 100 MHz resonator is smaller in size than the other resonator. A relatively 
high frequency flicker is therefore not surprising. 

The examples shown above indicate that, under the assumption of Sec- 
tions the oscillator frequency flickering is chiefly due to the fluctuation 
of the resonator natural frequency. 
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